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Parity-violating 3-jet observables for massive quarks to order oc^ in e + e 
annihilation 
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In this talk we discuss the calculation of the QCD corrections to parity-violating 3-jet observables in e + e~ collisions, keeping 
the full quark mass dependence M\. 



1. Introduction 

The measurement of the forward-backward asym- 
metry A FB for b-quark production in e + e~ collisions 
at the Z peak provides one of the most precise deter- 
minations of the weak mixing angle. It can be deter- 
mined with an error at the per mille level if the com- 
putations of A FB are at least as accurate as the present 
experimental precision of about two percent [|p. To 
achieve this precision one has to go beyond leading 
order in perturbation theory. Furthermore given the 
size of mass effects due to the non-zero fo-quark mass, 
it is also necessary to extend existing analyses to ac- 
count for non-vanishing quark masses. This is manda- 
tory if one studies top-quarks instead of b-quarks. Re- 
stricting ourselves to perturbative QCD the present 
knowledge is as follows. The order a s contributions 
had been computed first for massless [pf] and then for 
massive quarks ftpD- These calculations, which used 
the quark direction for defining the asymmetry were 
later modified [0-0] to allow the use of the thrust axis 
to define the asymmetry. The next-to-next-to-leading 
order (NNLO) coefficient was computed first in ref. 
[B], for massless quarks using the quark axis defini- 
tion. This result was recently corrected by a com- 
pletely analytical [|Kj] and by a numerical calculation 
rfTTI]. In ref. [ITTll the NNLO corrections were also 
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determined for the thrust axis definition. To complete 
the analyses of Afb at NNLO the corrections for mas- 
sive quarks are needed. To obtain these corrections it 
is convenient to compute the individual contributions 
of the parton jets [g], in particular to calculate sep- 
arately the contributions from 2-, 3-, and 4-jet final 
states involving a heavy quark Q. For the 4-jet contri- 
bution to Afb one needs only the known Born matrix 
elements for e + e~ — > 4 partons involving at least one 
QQ pair (see, e.g. ref. [|l2[|). The phase-space integra- 
tion of these matrix elements in order to get (ApB)4-jet 
for a given jet algorithm can be done numerically. 

Compared to (ApB^-jet the computation of 
(ApB^jet is much more difficult. Here the NLO 
corrections to the partonic subprocess e + e~ — ► QQg 
are needed. This contribution must be combined with 
the 3-jet contribution from subprocesses with four 
partons in the final state. Only by combining the two 
contributions a result which is free of soft- and mass 
singularities is obtained. In this talk we will present 
recent results [pi on the missing ingredients that are 
needed to calculate (ApB^-jet- 

The computation of (ApB^-jet for heavy quarks is 
beyond the scope of this work and remains to be done 
in the future. 

To fix our notation we discuss briefly in section 
the kinematics and the leading order results. In sec- 
tion we comment on the calculation of the one-loop 
corrections. In the following section some checks 
are described. We close with a conclusion in section 

B 



2. Kinematics and leading order results 

The leading order differential cross section for 

e + (p + )e-(p-) -> Q(k Q )Q(kQ)g(k g ) (1) 

can be written in the following formQ 

do 3 1 I ., , „. 

= --, — ^OuS (l+cos 2 r>)Fi 

dtydcostidxdx 4 (47i) 3 ' y ' 

+ (1 — 3cos f>) Fi + cos $7-3 + sin 2"r)cos <|)F4 



sin ■& cos 2<|)F5 + shrr}cos(])7'6 



with 



O pt = a(e + e -» f -» p+p ) = 



3s 



(2) 



(3) 



In Eq. (0), $ denotes the angle between the direction 
p_ of the incoming electron e~ and the direction kg 
of the heavy quark Q. The angle (]) is the oriented an- 
gle between the plane defined by e~ and Q and the 
plane defined by the quark anti-quark pair. The func- 
tions Fi depend only on the scaled cm. energies 



2k -k 



Q 



and x = 



2k ■ k? 



and on the scaled mass square 

»2 



m 
s 



(4) 



(5) 



of the heavy quark and anti-quark, where k = p + +p_ 
and s = k 2 . Note that the leading order differential 
cross section for a final state with n > 4 partons can be 
written in an analogous way in terms of the two angles 
"&, (j) and 3n — 7 variables that involve only the final 
state momenta. In higher orders absorptive parts give 
rise to three additional functions 7-7,8,9 (cf. ref. [|l3|]). 
The four functions F 1 2,4,5 are parity even. In partic- 
ular, the 3-jet production rate is determined by F\. In 
the case of massive quarks next-to-leading order re- 
sults for F\ were given in references [|l4|— pi|] . The two 
functions Fj and F(, are induced by the interference of 
a vector and an axial-vector current. In particular, us- 
ing the quark-axis, the 3-jet forward-backward asym- 
metry is related to F3 in leading order in the following 

4 We neglect the lepton masses and do not consider transversely 
polarized beams. 



way: 

(^Fb) 3-jet 



J da QcutQcut) sgn(cos(fl)) 

J da Q cut (y cut ) 
3 Jdxdx 8 cut (y cut ) F-j(x,x) 
8 Jdxdx Ocut^cut) F\ ( x ,x) ' 



(6) 



where 8cut(ycut) defines, for a given jet finding algo- 
rithm and jet resolution parameter y cat , a region in the 
(x,x) plane. The electroweak couplings appearing in 
Fjfi can be factored out as follows: 

^3,6 = 

-g$(l-l + l-)[Q Q ReX(g e a -f(K,l-)g e v ) 

+g^\i\Hf^ + ^-)(gf+gf)~-2 g e y a )\F X6 .{i) 

In Eq. (@), 



7/ -22/ sin 2 9 W , 



= T, 



3 ' 



1 
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4 sin 2 Qw cos 2 Qw s ~ m z + i m z^z ' 

/ — — A.4- 



1-A--A-4 



(8) 



where 7V is the third component of the weak isospin 
of the fermion /, Q\y is the weak mixing angle, and A, T 
denotes the longitudinal polarization of the electron 
(positron). In next-to-leading order in 0^, additional 
contributions to F 3 (, with electroweak couplings dif- 
ferent from those in Eq. (R) are induced which we will 
not consider here. They are either proportional to Im% 
and thus suppressed formally in the electroweak cou- 
pling or generated by the triangle fermion loop dia- 
grams. The contribution of the triangle fermion loop 
is gauge independent and UV and IR finite. It was 
calculated some time ago in ref. p2[]. 

The functions F^^(x,x) may be expressed in terms 
of functions h(,, hi which appear in the decomposi- 
tion of the so-called hadronic tensor as performed for 
example in references J23- 25 1 . 



F3 = ~[\/x 2 -4zh( l (x,x) 
+ \J x 1 -4zcos-r> e g/i7(jc, 



F 6 



= — vi 2 — 4zsin'd i 



QQ 



hi(x,x), 



(9) 



where ■&qq is the angle between Q and Q in the cm. 
frame and we have 



2(1 — x — x + 2z) + xx 
^ Vx 2 - 4zvx 2 - 4z 



(10) 



It would seem pointless to trade F^^ for h^j if it were 
not for the relation (which follows from CP invari- 
ance): 



hq{x,x) = —he(x,x). 



(11) 



To calculate the functions F3.6 it is thus sufficient to 
determine the function h(,{x,x). In the one-loop cor- 
rections to h(, one encounters both ultraviolet (UV) 
and infrared (IR) singularities. Regulating the UV 
as well as the IR singularities by continuation to d = 
4 — 2e space-time dimensions, we need the Born re- 
sult in d dimensions: 



/?6 {x,x) = \6%a s {N-l)B 



x 2x - (xx + x^ + 2-4x)e-4' 



with 
B 



1 



(12) 



(13) 



and the scaled gluon energy 

Xo— -=2-x-x. (14) 

s 

Note that the terms proportional to e in Eq. ( JI4 ) de- 
pend on the prescription used to treat 75 in d dimen- 
sions. To derive the above equation we have used the 
prescription 



W5^(l-^C F )^p Y5T PyV [26|?|. (15) 



3. Virtual corrections 

Given the vast knowledge on how to perform one- 
loop calculations the derivation of the one-loop cor- 
rections to h(, is in principle straightforward. So we 
restrict our discussion of the virtual corrections to 
some technical remarks and the final results. 

We used the background field gauge [28 2^] with 
the gauge parameter set to one. In this gauge the 
three-gluon vertex is simplified which leads to a re- 
duction of the number of terms encountered in in- 
termediate steps of the calculation. Furthermore we 



used the Passarino-Veltman reduction procedure 
to reduce the one-loop tensor integrals to scalar one- 
loop integrals. As mentioned earlier we regular- 
ized both UV and IR singularities by continuation 
to d = 4 — 2e space-time dimensions. Note that the 
't Hooft-Veltman prescription to treat 75 in d dimen- 
sions give rise to an additional 'finite renormalization' 
(see Eq. (|3])) to restore the chiral Ward identities. 

The finite contributions of the loop-diagrams are 
too long to be reproduced here. They are given explic- 
itly in ref. [[J]]. Here we give only the results for the 
UV and IR divergences and the contribution from the 
renormalization procedure. In the following we de- 
note with m always the mass parameter renormalized 
in the on-shell scheme and with oc 4 the strong cou- 
pling in the usual MS scheme evaluated at the renor- 
malization scale fj. 

Keeping only the pole-part of the integrals we ob- 
tain the following result for the UV singularities: 



1 virt., UV div. 



rr- 



1 / m 2 



L F h 6 



6 ' 1 eV4V7 I -" 

f 24a 2 (N 2 - l)C F 8h 6 + 9,a 2 (N 2 - l)C F Be ( 1 +x 

-2x+(l-x)-+z 



+ (2-3x)(l-x))+ 8i 

(l-x) 2 
+ 4zB(2 [ —^+x g ) 
1 —x 



2B({x 2 -2x)(l-x) 
x 2 -4z 



x(l —x) 



1 6a;C>fie( — (xx - Ax + 2 + x z ) 



Xo 



(16) 



with 
8h 6 



2zx g B(x ~ x)) + 0{£ 2 ) 



— 2 (x —3x — xx+ 3xx + 2x 



4x + 2) - x g (x 2 +xx-4x+ 2)e 



4 (-5x-3x + 2x 2 + 4+x 2 +xx)z 2 I, (17) 

1 —x 



and 



'1 = - (1 (fl^P (*(*+*) - 2 (! -*«))• < 18 ) 



The usual one-loop factor rp is given by: 

r(l+e)r 2 (l-e) 



rr = 



T(l-2e) 



The IR divergent contributions from the loop- 
integrals are given by 



i virt.. IR div. 



a, /VV ,lo 
2S ^J ^ 
1 1 



---ln(%-Iln(% 
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e z e ' m 

2 

ln(co) 



JV e 1 - co 2 
with tu = 2k( ■ kj and 



(20) 



1- i/l 



co = 



4z 



jc-Hc-1 



1 



4z 



(21) 



x+x-1 



Finally, the entire contribution from renormalization 
(including the Lehmann-Symanzik-Zimmermann 
residue and the 'ys-counterterm') is given by: 



/ virt., ren. 
h 6 



2^' TCF V4^ 2 



2 \ - £ 1 

m \ 1 



h LO 
"6 



24at(N z ~l)C F r r 



2 \ - £ 1 

m \ 1 



4K^ 



-8/j 6 
e 



0C S / m 2 

'IT T 



47T \ 47t^ 2 



- <UC f 



11 



-n)--N\\h 



t,LO 




2 1 1 . . ,m . 2 n m 2 i , 

- B/ - T A0ln( ? ) + 3lln^-8C, K 



-32a? (N 2 - l)C F 8h 6 - —CphY 3 



LO 



exactly by the UV singularities from the renormaliza- 
tion procedure. This cancellation is thus an excellent 
H9) check of the coefficients of the UV divergent inte- 

grals, namely the one- and two-point integrals. We 
checked also that the finite renormalization due to the 
't Hooft- Veltman prescription restores the chiral Ward 
identities. 

According to the Kinoshita-Lee-Nauenberg theo- 
rem [ pl| , p2[ ] the remaining IR divergences must be 
cancelled by the real corrections. To check this can- 
cellation we have calculated the singular contribu- 
tions from real emission by using a modified version 
|[5[||] of the phase space slicing method QfBJ]. By 
this cancellation the coefficients of the IR divergent 
triangle integrals and the box integrals are checked. 
The explicit results for the singular contributions of 
the real corrections are given in ref. [jjj], 

Furthermore we have verified that we are able to 
reproduce the coefficients of the loop-integrals in the 
massless result J23|-|25| ] . We found agreement with 
Ecl (4. 10) of reftlfTln the case of Eq. (4.9) of ref. 
[|4j we found agreement only up to an overall factor 
(—1). This is most probably caused by a typographi- 
cal error which is also present in ref. [£3j], but not in 
ref. pa] with which we fully agree. 



5. Conclusions 

In this talk we have presented the necessary in- 
gredients for calculating virtual corrections to parity- 
violating three jet observables. Together with the ex- 
isting results on the real contribution [Q, in partic- 
ular the singular contributions m\, these results allow 
for the calculation of parity-violating three-jet observ- 
(22) ables to next-to-leading order accuracy. 



Note that the singularity in the third line of Eq. (E2h is 
a collinear singularity and not a UV singularity. 

4. Checks 

In this section we discuss some checks to assure the 
correctness of our results. Given the fact that the fi- 
nite results are given in the form (coefficients) x (one- 
loop integrals) and that the loop integrals are known, 
it is sufficient to check the coefficients. As one can 
see from the comparison of Eq. ( |lq ) and Eq. ( |22| ) the 
UV divergences from the loop diagrams are cancelled 
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